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ABSTRACT 

We propose a new approach to measure differential radial velocities, mainly for single-lined 
spectroscopic binaries. The proposed procedure - TIRAVEL (Template Independent RAdial 
VELocities) - does not rely on a prior theoretical or observed template, but instead looks 
for a set of relative Doppler shifts that simultaneously optimizes the alignment of all the 
observed spectra. We suggest a simple measure to quantify this overall alignment and use its 
maximum to measure the relative radial velocities. As a demonstration, we apply TIRAVEL 
to the observed spectra of three known spectroscopic binaries, and show that in two cases 
TIRAVEL performs as good as the commonly used approach, while in one case TIRAVEL 
yielded a better orbital solution. 



Key words: methods: data analysis - methods: statistical 
techniques: spectroscopic - binaries: spectroscopic 



techniques: radial velocities 



1 INTRODUCTION 



Since the seminal works of ISimkinl J 19741) and iTonrv & Davisl 
the cross-correlation technique to measure astronomical 
Doppler shifts has become extremely popular. The advent of dig- 
itized spectra and computers rendered it the preferred method 
in astronomical fields that require measurement of radial ve- 
locities (RVs) from observed spectra, ranging from binary and 
multiple stellar systems to cosmology. In recent years, improve- 
ments in the precision of RVs measured through cross-correlation 
opene d the way t o the detection of many extrasolar planets (e.g., 
iMavor & Ouelozl 19951) . 

The cross-correlation technique is conceptually simple and 
can be presented in an intuitive manner. It searches for the velocity 
shift of the observed stellar spectrum which maximizes the cor- 
relation with a predetermined template. The template is assumed 
to adequately represent the object spectrum with no Doppler shift. 
The properties of the cross-correlation technique have been stud- 
ied extensively, in an effort to improve its precision and overcome 
its few limitations, and various methods ha ve been suggested to 



estimate its precision (e.g. 



l (e. 

iMurdoch & Hearnshawll991 



Tonrv & Davisl 1 19791 IConne3ll985k 
iBouchv et all200ll) . 



The precision of the individual RVs depends on several fac- 
tors, mainly on the amount of spectral information and the signal- 
to-noise ratio (S/N) of the spectra. Another crucial factor is the 
quality of the template spectrum, and the extent to which it rep- 
resents the actual stellar spectrum. Common practi ce is to use a 
theor etically calculated spectrum as template (e.g.. lLatham et all 
2002), a high S/N exposure of anot her similar star or the observed 
star itself (e.g.. lHowarth et alJl997h . 



The theoretical synthetic spectrum might have some system- 
atic differences relative to the actual stellar spectrum. This happens 
not only because the theory of stellar atmospheres is still not perfect 
and the spectral line lists are not complete, but also because every 
star is somewhat different from other stars, by its distinct abun- 
dances of the various atoms and ions in particular. Similar problems 
occur when an observed spectrum is used as template. An observed 
spectrum, either of the same star or of a different star, introduces 
random noise to the template, thus alienating it again from the true 
stellar spectrum. 

We present here a new approach, TIRAVEL, to derive the 
relative velocity of stellar spectra without using any template. The 
approach is based on the simple idea that when two spectra of the 
same star are available, one can use cross-correlation to measure 
a relative Doppler shift between the two spectra, each spectrum 
effectively acting as a template for the other one. TIRAVEL is a 
generalization of this approach for more than two spectra. We will 
show that for large number of spectra, this approach is equivalent 
to having a high-S/N effective template. Section|2|introduces the 
approach, while Section|3|applies it to three real test cases. We dis- 
cuss the test results and the applicability of TIRAVEL in Section 

m 



2 TIRAVEL 

TIRAVEL is basically a multi-spectral generalization of the con- 
ventional cross-correlation technique. When we cross correlate two 
spectra we scan a set of potential shifts between the spectra and 
score them by the degree of mutual similarity between the two 
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shifted spectra. This score is based on the correlation coefficient, 
which has been shown to be an o ptimal similarity measure under 
various assumptions fe.g. JZuckej 120031) . 

Suppose we wish to estimate the relative Doppler shift be- 
tween two spectra fi(n) and f2(n). The two spectra are described 
as functions of the pixel number - n, where n = A In A + B. Thus, 
the Doppler shift results in a uniform linear shift of the spectrum 
iTonr v & DaviJ 19791) . For simplicity we assume the spectra were 
already 'continuum-subtracted' and normalized, i.e.: 



= 



(i) 
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The cross-correlation function is then defined as: 

Ri2(s) = ^^fi(n)f 2 (n- 



s). 



This expression is essentially the correlation coefficient between 
the two spectra after they have been shifted according to the trial 
relative shift s. A higher value of Ri%(s) corresponds to bet- 
ter agreement between the two sequences. Furthermore, it is well 
known that |i?i 2 (s)| < 1. 

TIRAVEL generalizes the cross-correlation procedure to the 
case of a set of K > 2 spectra by considering the correlation ma- 
trix Rij (si, S2, ■ ■ ■ , Sk)- The entry Rij is the cross-correlation be- 
tween the i-th and j-th spectrum, as a function of their relative shift 
Si — Sj. In order to measure the overall agreement for a given set of 
shifts we propose to use Am - the largest eigenvalue of the matrix 
Rij (si, S2, • ■ • , s K )- We will show that this estimate reduces to the 
correlation coefficient for K = 2 

The rationale behind this measure is simple and intuitive. At 
the correct alignment, all the spectra are supposed to be easily mod- 
eled by one principal component, or template. The highest eigen- 
value of the correlation matrix measures the degree to which we can 
describe the spectra by one single principal component, as is well 
known from Prin cipal Component Analysis (PCA) theory (e.g., 
lMurtagh&Heckll987l) . The eigenvector of Rij, corresponding to 
the eigenvalue Am, can be used to produce the principal component 
itself, or in other words, the effective template TIRAVEL 'used' 
in order to produce the RVs. The elements of this eigenvector, with 
proper normalization, are the optimal weights with which the ind- 
vidual spectra should be summed to produce this 'template'. This 
is especially important in cases where the S/N ratios are varying 
considerably among the observed spectra. 

The value of the pairwise correlation coefficient can be inter- 
preted easily and intuitively, since values closer to unity mean bet- 
ter agreement between the two sequences, and values close to zero 
mean poorer agreement. We suggest to transform Am to obtain a 
similar behaviour: 

Am - 1 

p= lc=T' 

TIRAVEL essentially seeks the maximum of p as a function 
of the K relative shifts (si, S2, • • • , Sk )■ Each trial set of the K 
relative shifts requires the calculation of the corresponding corre- 
lation matrix Rij . Therefore, an important preliminary step in the 
computation should be the calculation of the K(K — l)/2 cross- 
correlation functions, corresponding to all the pairs of different 
spectra. The matrix elements Rij(si — Sj) are then simply sam- 
pled from the pre-calculated cross-correlation functions. 

Note also that the elements of the correlation matrix Ra de- 



Figure 1. Results of applying TIRAVEL to sets of K simulated spectra 
with zero Doppler shift. The plots present the mean RMS of the velocities 
produced by TIRAVEL. 



pend only on pairwise differences of shifts, i.e., relative shifts. 
It is therefore insensitive to zero-point shifts which are common 
to all the observed spectra. In that respect, the shifts derived by 
TIRAVEL are only relative shifts, with effective K — 1 degrees of 
freedom instead of K. 

After optimizing over {s\,s% , s K ), which c an be done by 
any optimization procedure (e.g.. IPress et aljfl992l) . we obtained 
the final velocities by fitting parabolas to the peak, similarly to the 
common practice in regular cross-correlation. 

Since the new procedure yields also an effective 'template', 
produced by summing all the individual exposures at the right shifts 
with optimal weights, we expect the quality of this effective tem- 
plate to improve as the number of observed spectra increases. Bet- 
ter template should yield more accurate individual velocities. We 
therefore expect the precision of the TIRAVEL velocities to im- 
prove as a function of the number of spectra. To test this we ran 
TIRAVEL on different sets of simulated spectra and checked the 
resulting velocities. We generated K spectra by using a synthetic 
spectrum of a G star to which we added normal noise at a speci- 
fied S/N. No shift was implemented in producing the spectra, and 
therefore the scatter of the velocities derived by TIRAVEL can 
serve as a measure of the precision of the approach. We repeated 
the procedure 50 times for each S/N and K. 

Figure Q displays the results of these simulations. The plot 
shows the mean RMS of the velocities obtained by TIRAVEL as 
a function of the number of spectra analysed, for three different 
S I 'N -ratios. As expected, the RMS converges to a constant value, 
when the 'effective template' becomes similar to the true template. 
It seems that the convergence rate is not very sensitive to the S/N, 
and that when 10-15 spectra are available, the ultimate perfor- 
mance is already achieved (assuming they all share about the same 
S/N). This behaviour is easily understood, since once the template 
is significantly less noisier than the individual spectra, the main 
factor which determines the velocity precision is the S/N of the 
observed spectrum. This happens already when a moderate number 
of spectra are co-added to produce the template. 

The simulations also indicate the computational burden of 
TIRAVEL. Using a very naive optimization scheme, on average 
about 10 iterations were needed to reach the stopping criterion of 
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less than 0.01 km s _1 RMS change in the velocities, for K = 30. 
We wrote our code in MATL AB and ran it on a Pentium Pro Linux 
machine, where those 10 iterations, including the pre-calculation 
step, lasted about 20 seconds. 



3 REAL TEST CASES 

We chose to demonstrate the potential of the new approach on 
three single-lined spectroscopic bin aries (SB Is ) which were al- 
ready solved and published (Latham et al. 2002) as part of a sur- 
vey of 171 high-p roper-motion SB Is. The radial velocities in 
Latha m et alJ <2002l) were derived from s pectra obtained using the 
CfA Digital Speedometers ( Lath amll98a. fl992). In order to obtain 
the radial velocities Latham et al. used an extensive grid of spec- 
tra computed using the model atmospheres of the ATLAS 9 code, 
developed by Kurucz (Morse & Kurucz, in preparation). 

The spectra we used for the TIRAVEL analysis were the 
same observed spectra that Latham et al. used for their analysis. 
Using the velocities produced by TIRAVEL we have recalculated 
a least-square best-fitting orbital solution. The RMS of the residu- 
als of the new orbit serves as a criterion to compare with the previ- 
ously published orbit. Furthermore, the effective template derived 
by TIRAVEL indicates the amount of spectral information that can 
be extracted from the spectra without using a template. 

3.1 G72-59 

Latham et all <2002l) solved G72-59 with 32 spectra and derived a 
period of about 88 days with a small eccentricity. Figure|2|shows 
three of the 32 spectra, to demons trate their S/N . TableQlists both 
the orbital elements that appear in Latham et al. (2002) and the or- 
bital elements obtained from the TIRAVEL velocities. Note that 
the TIRAVEL analysis reveals no information about absolute ra- 
dial velocities, and the value of the center-of-mass velocity ob- 
tained in analysing TIRAVEL velocities is not very informative. 
Nevertheless we list it for completeness. 

TableQshows that the two solutions are consistent. The resid- 
ual RMS of the TIRAVEL solution is somewhat larger than the old 
one, and so are the errors of the elements. This is on e of the cases 
where the TIRAVEL solution is a little worse than lLatham et alJ 
1 2002) solution. The two orbital solutions can be visually compared 
in Figure|3| Indeed, the differences between the two orbits are com- 
pletely negligible. Figure|4]presents the new effective template that 
is built by TIRAVEL, together with the synthetic template used in 
the original analysis. 

The synthetic template shown in Figure |4| was shifted accord- 
ing to the difference between the two derived center-of-mass ve- 
locities, in order to match the derived TIRAVEL template. The 
similarity is striking, although the synthetic template evidently has 
a much higher S/N than the derived one. It is somewhat surprising 
and gratifying that the TIRAVEL results are almost as good as the 
results of the calculated spectrum, despite the fact that the effective 
spectrum of TIRAVEL looks so much noisier. 

3.2 G178-27 

According to lLatham et alJ 120021) . G178-27 has a low amplitude 
orbit (K — 3.4km s _1 ), with an eccentricity of about 0.4. Its pe- 
riod is about 81 days, and t he orbital solution was obtained using 
31 exposures. lLatham et alJ derived very low metallicity for G178- 
27, of [Fe/H] = —2.0. Figure [5] presents three sample spectra of 
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Figure 2. Sample spectra of G72-59. The times of the presented exposures 
are JD = 2447053.7466 (upper panel), JD = 2447163.5952 (middle 
panel) and JD = 2449678.7457 (lower panel). 



Table 1. Comparing the two orbital solutions for G72-59 
Parameter Latham et al. 2002 TIRAVEL 



p 


[days] 


87.754 ±0.026 


87.755 ±0.031 


T 


[JD] 


2447 701.7±1.4 


2 447 701.8 ±1.5 


e 




0.168±0.017 


0.169 ±0.019 


Ul 


[°] 


260.8 ±5.4 


261.1±6.1 


K 


[kms" 1 ] 


15.33 ±0.24 


15.30 ±0.28 


V 


[kms" 1 ] 


11.78±0.18 


5.26 ±0.20 


&o-c 


[kms^ 1 ] 


0.94 


0.96 



G178-27. One can see the relative scarcity of the spectral informa- 
tion compared with G72-59, due to the very low metallicity. 

The two orbital solutions are listed in Table and shown on 
Figure [6] Once again, the new solution is very similar and consis- 
tent with the previous one, while its residuals are slightly smaller. 
The good choice of synthetic template bv lLatham et alJ is clearly 
seen once the TIRAVEL effective template is derived, as shown in 
Figure^ 

3.3 G48-54 

The star G48-54 is the best case considered here to demonstrate the 
potential of TIRAVEL. The previous solution by Latham et al. has 
a period of 22.6 days and an eccentricity of 0.3, obtained with 30 
exposures. The star G48-54 is the star with the latest spectral type 

Table 2. Comparing the two orbital solutions for G 178-27 



Parameter Latham et al. 2002 TIRAVEL 



P [days] 81.18 ±0.37 81.18 ±0.35 

T [JD] 2 447 436.5 ±1.7 2 447436.1 ± 1.8 

e 0.431 ±0.046 0.428 ±0.045 

w [°] 264.5 ±8.5 258.8 ±9.1 

K [kms" 1 ] 3.42±0.18 3.44±0.18 

Vo [kms- 1 ] -180.30 ±0.12 -1.11 ±0.13 



oo-c [kms" 1 ] 0.61 0.56 
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Figure 3. Upper panel: The orbital solution of G 72-59 using TIRAV EL 
velocities. Lower panel: The old orbital solution by Lath am et alJ 
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Figure 4. Upper panel: The effective template spectrum obtained by 
TIRAVEL for the spectra of G72-59. L ower panel: the Kurucz spectrum 
used as template by Latham et al. 1 2002), shifted to match the TIRAVEL 
effective template. 

presented in lLatham et al with an estimated temperature of 

3750 K. Figure[8]presents three sample spectra of G48-54. 

As can be seen in Table [3] the TIRAVEL solution is better 
than the old published solution. Because of the large RV ampli- 
tude, the differences between the two solutions can hardly be seen 
in the phase-folded RV curve (Figure We therefore present in 
Figure [TOl the residuals after subtracting the orbital solution in the 
two cases, on the same scale. The improvement is easily seen. Fig- 
ure \H\ suggests the reason for the improvement. The spectrum of 
this late-type star is more complex than most of the other stars in 
the sample, because of its lower temperature. We suggest that the 
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Figure 5. Sample spectra of G178-27. The times of the presented exposures 
are JD = 2446928.9137 (upper panel), JD = 2447494.0177 (middle 
panel) and JD = 2447702.7553 (lower panel). 
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Phase 

Figure 6. Upper panel: The orbital solution of G178-27 using TIRAVEL 
velocit ies. Lower panel: The previous orbital solution by Lath am et alJ 
120021) . 

spectral modeling of this star is not complete, probably due to the 
notorious difficulties associated with modeling later-type stars. 

Close examination of Figures l9l and 1 1 01 suggests that a single 
measurement, at a phase of about 0.3, may be 'pulled' by the sec- 
ondary spectrum, at a phase where the primary and secondary ve- 
locities are very close. This outlying velocity stands out even better 
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Figure 7. Upper panel: The effective template spectrum obtained by 
TIRAVEL for the s pectra of G178-27. L ower panel: the Kurucz spectrum 
used as template by Latham et al. 1 2002), shifted to match the TIRAVEL 
effective template 
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Figure 8. Sample spectra of G48-54. The times of the presented exposures 
are JD = 2445770.8684 (upper panel), JD = 2446861.8126 (middle 
panel) and JD = 2446959.6461 (lower panel). 



in the TIRAVEL velocities, because of the smaller RMS. Maybe a 
somewhat better procedure would exclude this point from the anal- 
ysis. 



4 DISCUSSION 

We have presented here a new technique to analyse the spectra of 
SB1 systems. TIRAVEL is independent of any assumption about 
the primary spectrum, which is not included in the actual observed 



Table 3. Comparing the two orbital solutions for G48-54 



> 

or 
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or 




Parameter 




Latham et al. 2002 


TIRAVEL 


P 


[days] 


22.6057 ±0.0049 


22.6020 ±0.0042 


T 


[JD] 


2 446 729.39 ±0.16 


2446 729.19±0.15 


e 




0.287 ±0.012 


0.257 ±0.010 


Ul 


[°] 


320.4 ±2.8 


316.2 ±2.6 


K 


[kms -1 ] 


25.74 ±0.32 


25.44 ±0.27 


Vo 


[kms- 1 ] 


0.84 ±0.24 


9.55 ±0.20 


&o-c 


[kms" 1 ] 


1.20 


0.89 



Figure 9. Upper panel: The orbital solution of G48-54 using TIRAVEL ve- 
locities. Lower panel: The previous orbital solution by Latham et al. l2002tl . 



data. Instead, TIRAVEL makes use of the observed spectra to de- 
rive the radial velocities of the primary. This approach assumes that 
there is only one component in the observed spectra and that they 
are free of atmospherical or mechanical systematic features. Such 
features might reduce the efficiency of the technique, in particular 
compared to the use of theoretically calculated templates. 

We have applied TIRAVEL to three SB1 cases, where a satis- 
factory solution had already obtained with theoretically calculated 
templates. TIRAVEL improved significantly the quality of one of 
the three solutions, in terms of the precision of the orbital elements 
and the RMS residuals. Obviously, the advantages of TIRAVEL 
will be most pronounced when there is no good theoretically cal- 
culated template available. In the presented cases, the template grid 
that lLatham et alJ 120021) used was very extensive, and we could 
not expect to obtain very large improvements, except at the low- 
temperature frontier of the grid. One immediate implication of 
TIRAVEL is the analysis of the spectra of cataclysmic variables 
or some T Tauri stars, where the high-temperature accretion disc 
has an emission spectrum which in not easy to model. 

The version of TIRAVEL we presented here is specifi- 
cally tailored for SB Is, where there is only one RV to be mea- 
sured in each exposure, resulting in a simple shift of the spec- 
trum. In the case of double-lined spectroscopic binaries (SB2), 
a more elaborate scheme is needed, one that has to sepa- 
rate the effects of the orbit on the two component spectra. 
Such a sc heme may make use of ideas related to Doppler to- 
mography ( Bagnuolo et al. 1992) or two-dimensional correlation 
(TODCOR lZucker&MazehUl994h . One possible generalization 
of TIRAVEL may simply apply it iteratively - first solving for 
the primary component, subtracting the resulting template from the 
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Figure 10. Upper panel: The residual RVs after subtracting the best-fitting 
orbital solution from the TIRAVEL velocities of G48-54. Lower panel: 
The residual RVs after subtracting th e previous best- fitting orbital solution 
from the velocities published by Latham et al. 1 2002). 
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Figure 11. Upper panel: The effective template spectrum obtained by 
TIRAVEL for the spect ra of G48-54. L ower panel: the Kurucz spectrum 
used as template bv lLatham et alJ i2002h . shifted to match the TIRAVEL 
effective template 



spectra, then applying TIRAVEL to the subtracted spectra again 
to retrieve the information related to the secondary component, and 
repeat the process iteratively. Similar procedures were applied by 
iMarchenko et alHl998h and lGonzalez & Levatd <2006t) . 

Note th at spectral disentangling codes, like KOREL jH adrava 
Il995lll997h can be applied also to SBls fe.g.. lSaad et alJl2005h 
However, they usually impose an orbital model, treating any de- 
viation from the model as noise. TIRAVEL does not incorporate 
any orbital model (or any other source of RV variation), and pro- 
duces the best relative RV estimates one can get without a template. 
Therefore, the TIRAVEL approach enables us to detect any addi- 
tional source for RV variation, a third companion in particular. 

In order to be applicable to spectra obtained by the modern 
echelle spectrographs and other multi-order spectra, TIRAVEL 
has to be generalized, probably along lines similar to the ones that 
were used to generalize TODCOR to multi-order spectra iZuckeJ 
120031) . Thi s modification was cruc ial in detecting the planet in 
HD 41004 feucker et alj2003ll2004h . 

In principle, improvements in RV precision, like the one of- 



fered by TIRAVEL, may lead to the detection of small finite ec- 
centricities or additional stellar components. One area where the 
precision of radial velocities and orbital solution is of utmost im- 
portance is the search for extrasolar planets, where the orbital RV 
amplitude is usually extremely small, of the order of tens of me- 
ters per second. Effective improvements in the precision of radial 
velocities may lead to an improved planet-detection capability. In 
general, the improved precision also means a better utilization of 
the instrumental potential and thus a more efficient observing time 
usage. 
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